Quantum to classical crossover in the 2D easy-plane XXZ model 
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Ground-state and thermodynamical properties of the spin- 1/2 two-dimensional easy-plane XXZ 
model are investigated by both a Green's-function approach and by Lanczos diagonalizations on 
lattices with up to 36 sites. We calculate the spatial and temperature dependences of various 
spin correlation functions, as well as the wave- vector dependence of the spin susceptibility for all 
anisotropy parameters A. In the easy-plane ferromagnetic region (—1 < A < 0), the longitudinal 
correlators of spins at distance r change sign at a finite temperature Jb(A, r). This transition, 
observed in the 2D case for the first time, can be interpreted as a quantum to classical crossover. 



I. INTRODUCTION 

The magnetic properties of low-dimensional quantum 
spin systems with spin anisotropy, such as the quasi-one- 
dimensional (ID) cuprates [1] and the quasi-2D high-T c 
parent compounds [2], are of growing interest. The S = 
1/2 XXZ model 
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denote nearest- neighbor (NN) sites; throughout 
we set J — 1) usually serves as the generic model for 
those systems. Recently, in the ferromagnetic (FM) re- 
gion (— 1 < A < 0) of the ID model a quantum-classical 
crossover in the longitudinal spin correlators was found 
by means of exact diagonalization (ED) [3] and a Green's- 
function theory [4]. For the XXZ model on a square 
lattice, an analytical approach to the spin susceptibil- 
ity taking into account the magnetic short-range order 
(SRO) at arbitrary temperatures does not yet exist. 

In this contribution the spin correlations in the easy- 
plane region — 1 < A < 1 of the 2D XXZ model are 
examined by both a Green's-function theory outlined in 
the Appendix and by exact finite-cluster diagonalizations 
of the model (1) on lattices with up to 36 spins using peri- 
odic boundary conditions. We mainly focus on the char- 
acteristics of a possible quantum to classical crossover 
in the FM regime. Moreover, for the first time, the com- 
plete wave- vector, temperature and A dependences of the 
static transverse and longitudinal spin susceptibilities are 
calculated. 



II. GROUND-STATE PROPERTIES 

In Fig. 1 our results for the magnetization m(A) are 
compared with available quantum Monte Carlo (QMC) 
data [5], where the ED/QMC data for the ground-state 
energy per site e(A) (inset) is taken as input for the 



Green's-function approach (Cfp = ^de/dA, C^ 
e/2-ACg). 
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FIG. 1.: Magnetization m and ground-state energy e of 
the 2D easy-plane XXZ model. 

As can be seen from Fig. 2, the short-ranged corre- 
lations calculated analytically are in excellent agreement 
with our ED data. Let us stress that the finite-size depen- 
dence of the ED data is almost negligible by going from 
a 32- to a 36-site lattice. At A = 1 the rotational sym- 
metry Cjj = 2C ZZ is visible. At the quantum critical 
point A = -1 we have C+T = 2C ZZ = 1/6 (cf. Eq. (14)). 

The non-analytical limiting behavior lirriA— >-i+ C zz = 
results from both the QMC [5] and ED data (obtained in 
the subspace with total spin projection S z = 0). 

The static spin susceptibilities Xq(A) are depicted in 
Fig. 3. In the FM region, for sufficiently low A values, 
Xq shows a maximum at q = being a precursor of 
the FM instability (in the zz-correlators) at A = — 1. 
Note that (Xq~) _1 — 0, reflecting the transverse long- 
range order (LRO) at T = 0, by Eq. (13) corresponds 
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0. In the antifcrromagnetic (AFM) region 



< A < 1 the maximum in x^ z at q = Q is indicative 
of the longitudinal AFM LRO at A > 1. 

Finally, in Fig. 4 we show the longitudinal spin-wave 
spectrum iv^ z (cf. Eq. (7)). For q = |q| <C 1 we have 
uj^ z = c z a z q, where the spin-wave velocity c zz increases 
with A over the whole easy-plane region. The minimum 
in at q = Q in the AFM region corresponds to the 
maximum in x* 2 (cf. Fig. 3 a) and reflects the increase 
of the longitudinal AFM SRO with A (see also C ZZ (A) 
in Fig. 2). 
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FIG. 2.: Transverse and longitudinal spin correlation 
functions at T — 0. Symbols denote ED results ob- 
tained for a 6x6 lattice. 



III. FINITE-TEMPERATURE RESULTS 

The temperature dependence of the short-ranged lon- 
gitudinal spin correlations is displayed in Fig. 5. Again 
the analytical results agree remarkably well with the ED 
data. In the FM region, for the first time in the 2D model, 
we observe the so-called "sign-changing" effect which was 
found numerically [3] in the ID model and later on repro- 
duced by our Green's-function calculations [4]. That is, 
at fixed separation r and with increasing temperature or 
at fixed temperature and with increasing r, C zz changes 
sign from negative to positive values. The temperature 
T (A, r) where C ZZ (T (A, r), A) = are given in Table I. 
As in the ID case, T at fixed A decreases with increasing 
r. However, compared to the ID case [4], our analytical 
results are in much better agreement with the ED data. 
The sign change of C zz may be interpreted as a quantum 
to classical crossover [3] because with increasing temper- 
ature the system behaves more classically, i.e., it becomes 
dominated by the potential energy (negative A term of 
the Hamiltonian favoring the parallel alignment of two 
spins). In the AFM region we obtain the expected al- 



FIG. 3.: Wave-vector dependence of the longitudinal (a) 
and transverse (b) static susceptibilities Xq at T = 0. 
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FIG. 4.: Longitudinal spin- wave dispersion ui^ along the 
major symmetry directions of the 2D Brillouin zone. 



ternating signs of C zz corresponding to the longitudinal 
AFM SRO. 

In Fig. 6 various susceptibilities a * q = 0, Q are 
plotted as functions of T and compared with numeri- 
cal data. For A = 0.5 the longitudinal and transverse 



2 



uniform susceptibilities are in reasonable agreement with 
the QMC results [5] and our ED data (the up- and down- 
turn at lower temperatures is a finite-size effect). The in- 
crease of x${T), the maximum near the exchange energy 
(J = 1), and the crossover to the Curie- Weiss law are 
due to the decrease of AFM SRO with increasing tem- 
perature. On the other hand, the staggered susceptibility 
Xq is enhanced as compared with Xq z by the longitudi- 
nal AFM SRO. In the FM region (Fig. 6 b, A = -0.5) 
the maximum in Xq z S where the analytical and numerical 
results yield nearly the same position, may be explained 
as a combined SRO and sign changing effect as discussed 
for the ID model in Ref. [4]. Contrary to the AFM re- 
gion, Xq is suppressed as compared with Xo z which is 
caused by the FM correlations above To. The tempera- 
ture dependence of Xo~ = may be explained again 
as a SRO effect. Here, the transverse FM SRO results in 
a spin stiffness against the orientation of the transverse 



spin components along a staggered field perpendicular to 

the z-direction, so that xt^c IS suppressed at low teln- 
et, ri 

peratures and exhibits a maximum. 

TABLE I. Temperature To(A;r) of the sign change in 
the longitudinal correlation functions C^ Z (T; A). The corre- 
sponding results obtained from ED of a 4x4 lattice are given 
in parenthesis. 
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FIG. 5.: Temperature dependence of the NN (a) and next 
NN (b) longitudinal spin correlation functions C* z ■ Sym- 
bols denote ED results obtained for a 4x4 lattice. 



FIG. 6.: Longitudinal and transverse static spin suscep- 
tibilities Xq as functions of temperature T for the 2D 
AFM (a) and FM (b) easy-plane XXZ models. 
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IV. SUMMARY 

To resume, we presented a Green's-function theory of 
magnetic LRO and SRO in the 2D easy-plane XXZ model 
which allows the complete calculation of all static mag- 
netic properties in excellent agreement with numerical 
diagonalization data. In particular, in the FM region we 
found a quantum to classical crossover in the longitudi- 
nal spin correlations. We conclude that our approach is 
promising to be applied to other anisotropic spin mod- 
els, such as the quasi-2D XXZ model for the parent com- 
pounds of high-T c superconductors. 



APPENDIX: GREEN'S-FUNCTION THEORY 

The spin susceptibilities x+~(u) = — ((S+; SZq)) u an d 
Xq Z (w) = — ((Sq, Si q )} w , expressed in terms of two-time 
retarded commutator Green's functions, are determined 
by the projection method, developed, for the XXZ chain, 
in Ref. [4]. Taking the two-operator basis (Si", *S q ) T and 



(S*, iS*) we obtain 



q 



4 LO Z — (luV) z 



; v = H — , zz, 



(2) 



with 



M+~ = -4[C+-(1 - A 7q ) + 2Cf *(A - 7q)] , (3) 



^q Z = -4C+-(l- 7 q) 



(4) 



Cl m = C», C+~ = (S+S-), C« = (S Z S Z ), r = ne x + 
me y , and 7 q = (cosq x + cosq y )/2. The spin correlators 
are obtained from Eq. (2) as 

1 . M v 



(5) 



where p(ojq) = (e w «' — The spectra w q , calculated 

in the approximations — S q = (w q ~) 2 S q and — S q = 

(w q z ) 2 S q introducing vertex parameters a" (i = 1,2), 
are given by 

«-) 2 = [(1 + 2a+-(C+- + 2C+-)](1 - A 7q ) 
+A(1 + 4a+-(C zz + 2C£)](A - 7q ) 



+2 a +-[C+-(4A 7 2 -A-3 7q ) 
+2C^(4 7 2 -l-3A 7q )], 

«) 2 = 2(1 - 7q )[l + 2a? (C+- + 2C+-) 
-2AafC+-(l + 4 7q )]. 



(6) 



(7) 



In the easy-plane region — 1 < A < 1, the long-range 
order at T = is reflected in our theory by = 

[Q = (tt,7t)]. Accordingly, the condensation part Ce jCJr 



is separated from C+ (cf. Eq. (5), and the magnetiza- 
tion to is calculated as 



TO 



= lyc+-e-^ r = c. 



(8) 



The parameters a\ (T) are determined from the sum rules 
C£ ~ = 1/2 and C§§ = 1/4. To obtain o£(T) we adjust 
Cio(T = 0) taken from our ED data and assume, as 
additional conditions for the calculation of X q z ( w ) an( l 
X q ~(w), temperature independent ratios 



R zz = 



a Z 2 Z (T) - 1 
af (T) - 1 



and 



R +- = a t~{T)-l for A>0 
> a+-(T)-l 



«2 + "(T)-l 

a\ z {T) - 1 



for A < , 



(9) 

(10) 
(11) 



respectively. For the discussion it is useful to perform 
the unitary transformation which rotates the spins on 
the sublattice B around the z-axis by the angle 7r, S, = 
U+S t U with U = UieB 2S i- We g et <S?' y = e iC * Ti S*' y , 
S? = Sf and 



(12) 



{id) 



Due to {A)h = for any operator A, we obtain the 

relations 



xtn{ UJ ) = xtA^)\ k = q-Q, 



C, 



r,H 



(13) 
(14) 

Xq/HM = Xq^(w), and C ZZ H = C zz ^. As shown in 
Ref. [4], the rotational symmetry at A = ±1 is preserved 
by our theory. 



[1] D. C. Johnston, in K. H. J. Buschow, editor, Handbook 
of Magnetic Materials, volume 10, Elsevier Science, (Am- 
sterdam 1997). 

[2] M. Matsumura, F. Raffa, and D. Brinkmann, Phys. Rev. 

B 60, 6285 (1999). 
[3] K. Fabricius and B. M. McCoy, Phys. Rev. B 59, 381 

(1999). 

[4] C. Schindelin, H. Fehske, H. Biittner, and D. Ihle, cond- 
mat/9911447. 

[5] Y. Okabe and M. Kikuchi, J. Phys. Soc. Japan 57, 4351 
(1999). 



4 



